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We present a methodology to renormalize the nucleon-nucleon interaction in momentum space, 
using a recursive multiple subtraction approach that prescinds from a cutoff regularization, to con- 
struct the kernel of the scattering equation. The subtracted scattering equation is solved with 
the next-leading-order (NLO) and next-to-next-leading-order (NNLO) interactions. The results are 
presented for all partial waves up to j = 2, fitted to low-energy experimental data. In this renor- 
malization group invariant approach, the subtraction energy emerges as a renormalization scale and 
the momentum associated with it comes to be about the QCD scale (Aqcd), irrespectively to the 
partial wave. 
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I. INTRODUCTION 

The nucleon-nucleon (NN) interaction in leading order corresponds to the one-pion-exchange potential (OPEP) plus 
a Dirac-delta, when considering an effective field theory (EFT) of nuclear forces based on a chiral expansion of the 
effective Lagrangian. This procedure was suggested by Weinberg [1] , with a recipe to infer the values of the strength 
of the Dirac-delta interaction in the ^So and '^S'l channels from the singlet and triplet scattering lengths respectively. 
Therefore, the singlet and triplet scattering lengths allows to fit the renormalized strengths of the contact interactions. 
This effective potential should be valid for momenta well below some typical momentum scale considered in quantum 
chromodynamics (QCD), such as the rho meson mass {rrip ~ 4 fm~^) [I], which implies in a cutoff at the momentum 
scale of this order or below, for the intermediate virtual propagation of the NN system. 

According to Weinberg's recipe, one should use naive dimensional analysis (NDA) to order terms in the potential, 
truncate it at a certain order, and then solve the Lippmann-Schwinger (LS) equation exactly. However, it has been 
shown that this recipe is not consistent with renormalization-group (RG) invariance, already at leading order [2-4]. 
The validity of Weinberg's power counting is questioned [2H1] in particular for the waves like ^Pq, where the singular 
and attractive pion tensor force is requiring more than two parameters even at LO. 

About a decade ago, an alternative way to renormalize the nucleon-nucleon interaction, for a singular potential was 
proposed in Ref. [5]. In an extension of |S], the approach was proved in Ref. [5] to be RG invariant, as it satisfies 
the corresponding non-relativistic Callan-Symanzik equations [?]■ In this approach, no cutoff is considered in the 
equations and/or interactions. Instead, it is introduced a subtraction point in the kernel of the Lippmann-Schwinger 
(LS) equation [8|, in order to reach a finite T— matrix. Our research group has considered singular contact interactions 
in the context of nuclear [51 [SI IHl HO] , atomic and general physics [TTJ[T2]. In such works, within specific renormalization 
procedures applied to few-body systems, scaling limits and correlations between low-energy observables emerge as a 
consequence of using singular contact (zero-range) interactions. 

As the approach was proved to be invariant under renormalization group transformations, after fitting the data, 
one can easily move the reference scale without affecting the physics. This flexibility in moving the reference scale is 
a big technical advantage of the present approach in relation to other ones. Indeed, the renormalization results for 
the observables should be independent of the specific scheme used. Our method was applied to the NN interaction 
with OPEP supplemented by contact interactions [H [SJ [TD] . 

The significant results obtained by several authors in the implementation of the EFT program for the two-nucleon 
system [T51[25] include a vast literature on the predictive power of the leading order term (OPEP plus delta) with a 
single renormalization momentum scale. Such calculation gave a sound basis for the renormalization program of EFT 
in the NN system. In particular, the OPEP background to the NN observables were carefully analyzed in Ref. [29j. 

Renormalization of the NN interaction in chiral effective theory in leading order (LO) , up to next-to-leading order 
(NLO), and to next-to-next-to-leading order (NNLO), in S, P, D waves, has been done with success in coupled and 
uncoupled waves [5[ \6\ ITO l l30H34] . The references [30j 131j apply the renormalization approach in coordinate space, 
while the works ^ [HI [IHl [Ml IM] use subtracted-resolvent two-body techniques in momentum space. The subtracted 
renormalization procedure considered in Refs. [331121] is essentially different from the one applied in Refs. [Sl[51[Tn] in 
the way to deal with the physical inputs introduced to heal the ultraviolet divergences. In short: in the renormalization 
strategy of [33] . one- folded subtracted equation and also a cutoff are used, implying in a window where the observables 
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are quite independent on the cutoff. The strategy of Refs. [SJinilin] is based on the ehmination of ultraviolet divergences 
and relies on multiple subtractions, without any additional cutoff parameter. 

The multiple subtraction technique that we are going to use in the present approach demands a strategy to construct 
the driving term of the corresponding scattering equation, which is generated by consecutive subtractions in the kernel 
at some defined energy scale. We denote a generic order of subtractions in the kernel by n, such that n — I stands for 
one subtraction. With the assumption that we have only 5— wave LO contact interaction, one subtraction is enough. 
In this case the subtraction scale is let to infinity and then driven towards a given finite value of the subtraction 
energy by solving a non-relativistic Callan-Symanzik (NRCS) equation [SUH]. The corresponding solution gives the 
LO driven term calculated at the reference subtraction energy where more contacts and TPE potentials are included. 
We introduce each order of the potential (LO, NLO, NNLO) in the driven term of the subtracted equation when the 
number of subtractions are enough to make it finite. For that, we use n = 3 in case of NLO, and n = 4 in NNLO. 
Unitarity is strictly kept. In this regard, we note that the lowest contact of the P— wave is introduced together with 
the TPEP NNLO. In the complexity of this subject, to have different approaches to tame highly divergent potentials 
inherent to Weinberg proposal, it is important to observe the coherence between the physical results obtained with 
higher order potentials and the corresponding contacts to learn subtle aspects of QCD that permeates the chiral 
effective expansion of the nuclear interaction. In other words, a trace of the QCD scale is expected to emerge in the 
fitting procedure for the nucleon-nucleon scattering data. And, as we are going to verify in our approach, this relic 
can be found in the scale associated with the subtraction point where the two-pion exchange potential at NLO and 
NNLO are introduced in the renormalization process. 

In the present work, going beyond the LO interaction and also considering P— and £)— partial waves, we show how 
to implement a multiple subtractive renormalization approach to obtain consistent results for the NN observables. 
We start by summarizing the methodology used to renormalize the nucleon-nucleon interaction using a recursive 
subtracted approach that was previously considered in Ref. j^jSKTU]. Next, we present the equations to be solved in 
order to obtain a finite T— matrix with NLO and NNLO interactions, taking into account partial waves up to j = 2. 
The results are fitted to low-energy experimental data. As it will be shown, the adjusted derivative contact terms 
dominates over the NLO and NNLO interactions in the S*— wave channels for momentum higher than about 1 fm"-'^. 
It is about the same, in the case of nonzero angular momentum. By including higher contact derivative terms, the 
potential is fitted to low energy phase-shifts up to about 200 MeV laboratory energies. 

As we have verified, the subtractive-renormalization technique considered in Refs. [5j [6] is shown to be reliable also 
for highly-singular potentials up to NNLO and j = 2. These potentials have also been discussed in Refs. [551 [31]. 
In order to show that, we present results for the nucleon-nucleon phase-shifts up to j = 2, going beyond previous 
calculations within the present technique. We stress that our method cannot be naively confused with a Born 
approximation which, in principle, is not reliable for singular interactions. Indeed, we show that the unitarity is 
strictly kept along all the subtractive procedure. Moreover, it was already demonstrated in Ref. [S] that our multiple 
subtractive approach is fully renormalization-group (RG) invariant. 

We also show how to implement the multiple subtractive renormalization method for P— waves with contact inter- 
actions, respecting the phase-shift threshold behavior of these waves. By considering that, we supply the details not 
explicitly discussed in Refs. [5, |3j. Actually, P— wave channels were already studied in Ref. |34j with one subtraction 
and a sharp momentum cutoff; and, in Ref. |31j (and references therein), in coordinate space, even considering higher 
partial waves. 

The number of recursive steps required to renormalize the interaction depends on how the potential diverges. For 
instance, the leading order requires only one step, the next-to-leading order requires three steps and the next-to- 
next-to-leading order requires four steps. The method has been extended for a generic number of recursive steps 
[6] and applied to one-pion exchange plus contact interactions with three steps |10) . Our present renormalization 
technique can be used to organize and implement calculations where one-pion exchange is treated non-perturbatively. 
Higher order terms from the chiral expansion of the nucleon-nucleon force can also be treated within this method. A 
further advantage of our method can be particularly appreciated in the ^Sq and ^Pq channels, in which we show the 
phase-shift behaviors for increasing cutoffs, exemplifying the cutoff independence as it goes to infinity. In conception, 
our method differs from the works presented in Refs. [33JI31] and [3D1I3T] which also dealt with simpler effective chiral 
expansion in the way suggested by Weinberg. In other approaches, such as the one given by Barford and Birse (35], 
the higher chiral order are evaluated in distorted-wave Born approximation. However, in this case the analysis of the 
renormalization-group is quite involved. In Ref. |19j . it was considered a systematic chiral procedure with a smooth 
cutoff applied to the potential in order to obtain the observables from the originally divergent scattering equations. It 
has also been suggested that subleading orders of the potential should not be fully iterated, but treated in (distorted- 
wave) perturbation theory (see e.g. |35j and |36j ) . In our approach, we can iterate to all orders the potential and still 
keep the renormalization group invariance. 

This work is organized as follows. In section II, we briefly describe the subtractive formalism applied to the 
Lippmann-Schwinger equation, with the renormalization group equations. The formalism is followed by an analytical 
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example where the RG equation is solved for the P— wave case with contact interactions, where two subtractions 
are required. In section III, the renormalization formalism is applied to the NNLO potential, by considering recur- 
sive renormalization processes. The results for the NN phase-shifts are presented and discussed in section IV. Our 
conclusions are given in section V. 



II. SUBTRACTED LIPPMANN-SCHWINGER EQUATION 

The scattering equation can be written with an arbitrary number of subtractions (n) in the kernel [Gj 110) , that is 
useful when the potential has ultraviolet divergences with attractive nature e.g. — 1/r™ with m > 2 (for m = 2 the 
critical strength for collapsing the bound state is 1/4), and/or the potential includes contact terms (Dirac delta and 
it's derivatives). By considering our units such that H and the nucleon mass tojv are equal 1, we have the energy E 
given hy E = k'^. Within such units, the regularized subtracted form of the LS equation, in operatorial form, can be 
written as 



where, by having Hq as the free operator, we have 

/ ,,2 , P \ " 

G(o)(S), 



(1) 



(2) 



(3) 



We should note that, n is supposed to be the necessary maximum number of subtractions, with — /i^ the energy scaling 
parameter, to render finite results for the T— matrix. As we can verify in the above equations, the dependence on the 
energy starts to appear in y*^") for n>2. The formal solution of the regularized Eq. (hj), for V^*^") = V^"'\E; — /i^), is 
given by 



(4) 



By considering the operatorial expression given by Eq. ([s]) with n subtraction, in explicit momentum-space notation, 
for an arbitrary single partial-wave (after integrated the angular parts and with the assumption that the interaction is 
symmetric and not angle dependent), and by taking {q,P', fc^; — /x^) as the matrix element of the angular momentum 
projected operator, with 



dpp^ \p){p\ 



we obtain the partial-wave projected equation, which is given by 



(5) 



dq q '{q,q\k )^-— — j——V^>{q,p;k-~^i) 



ifJ.^ + q'- 



(6) 



Therefore, as an example, we can evolve in subtractions the OPEP plus delta, which is recognized to be renormalizable 
by fixing only one S*— wave observable in the coupled triplet and singlet states (usually the scattering length is chosen 
as the input for fixing the renormalized strength of the contact). In this particular case, when only one subtraction 
is enough to produce a finite T-matrix, V^(^^(--/i^) is replaced by the corresponding T— matrix at the point E = — /i^, 
such that, by defining T{~fi^; —jJ^) = T{—fi'^), we can obtain the following subtracted equation: 
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For the Dirac delta-potential, the matrix element in momentum space of T{—iJ?, ~A*^) is the renormalized strength. 

The subtraction energy — /z^ in Eq. ([l]) can be moved without changing the resulting T— matrix if the driving term 
is evolved by solving the Callan-Symanzik equation [5J [71 [TU] . The renormalization group equation for the driving 
term is briefly sketched below. One example of renormalization of a P— wave case, with our method, is also discussed 
for illustration. On the calculations of shallow P— wave states, with EFT theory applied for halo nuclei, see Ref. [57] . 

In addition, we remark a general feature of the subtraction procedure applied to attractive and repulsive power-law 
potentials in configuration space (r~™ with m > 2). In the repulsive case the subtraction is not trivially required to 
render finite the T— matrix, as the scattering wave function is damped in the classically forbidden region. In contrast, 
in the attractive case the ultraviolet divergence is actually met dynamically and may collapse the state. That is 
the case that a regularization has to be performed at the short range. In our method we deal with this problem, 
by using enough subtractions to make the T— matrix finite. Attractive and repulsive terms that we are discussing 
now appears as the chiral expansion is performed, as also discussed thoroughly by Valderrama and Arriola in their 
renormalization scheme in configuration space for calculation of the NN phase-shifts [30] . Instead, in our approach, 
subtractions are introduced to include the contacts, at the expense of finding a convenient subtraction point with the 
meaning of a reference low-energy scale. As we will show, indeed the fitting of partial waves phase-shifts up to j = 2 
gives a reference subtraction energy between -50 and -100 MeV. At this scale the NLO and NNLO TPEP potentials 
are identified as contributions to the driving terms as we are going to explain in detail in Sect. III. 

Note that the reference subtraction point is, in principle, arbitrary, and can be moved by evolving the potential 
through the RG flow equation. Once the fitting to the data is done, the reference scale can be changed arbitrarily as 
long as the driving term (starting potential) is evolved through the nonrelativistic Callan-Symanzik flow equation [7]. 
In fact, this implies that the RG invariance cannot be separated from the form of the evolved interaction. Next, we 
briefly discuss the RG equations, followed by an analytical example. 



The nucleon-nucleon observables are invariant under the change of the arbitrary subtraction point, therefore one 
can start at any convenient energy scale — /i^. However, the form of the driving term and its coefficients, which define 
the scattering amplitude are tied to the prescription used to define the renormalized theory. The key point of the 
renormalization group method is to change this prescription without altering the predictions of the theory [38j . 

The invariance of the T— matrix under changes of renormalization prescriptions, i.e., 



which is derived from the invariance of T{E; —fJ.'^), given by Eq. Q, with respect to the renormalization parameter 
fi^. The demonstration of the above can easily be done by considering both expressions that appear in Eq.Q. 
Equation Q substantiate the invariance of the renormalized T— matrix under dislocation of the subtraction point. 
Then, we observe that there is a non-trivial dependence on the subtraction point appearing in the driving term of the 
subtracted scattering equation, although the physical results of the model are kept unchanged. Thus, from now on, 
we drop the explicit dependence on /i^ in the T-matrix by writing T{E; —fi^) = T(E). 

The solution of ([s]) implies in a complicated evolution of 1/'"^ as fi changes. Not only the strengths of the interactions 
would change, but also the form of the driving term. The ultraviolet behavior of the driving term is not changed by 
the evolution in fi. The evolution should not be truncated as fi is varied in order to keep the T— matrix invariant. 
At different /i the potential /i^; E) has a complicate form from the solution of NRCS equation. Similarly, the 

evolution of renormalization group equations as introduced by Bogner, Kuo and Schwenk (3^] for the NN scattering 
does not truncate on certain operators as the cutoff is varied to keep the observables unchanged. 



A. Renormalization Group Equations 
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B. Example: Subtracted P— wave equations with contact interactions 

We illustrate the use of our multiple subtraction renormalization method by discussing the case of a one-channel 
P— wave problem with a derivative contact interaction V{p,p') — Ai p p' ■ The P— wave Lippmann-Schwinger equation 
with a momentum cutoff A is easily solvable, with the corresponding T— matrix given by 



r(p,p';fc^ 



1 2 



dq q 



1 



jr.2 _ q2 ie 

The coefficient Ai can be fixed by the scattering volume a, which gives 

Ai dTT a 



By replacing Ai in Eq. (10), we have 



T{p,p'; 



■ p-p 



n 



'^'^^ k^-q^+ie 



(10) 



(11) 



(12) 



which shows that this procedure is not enough to eliminate the cutoff dependence of Eq. (10). With this simple 



example, we can see that a renormalization procedure using only one subtraction, for a P— wave problem, is expected 
to be cutoff dependent, requiring a second subtraction. 

The fit of only the scattering volume is equivalent in our method to use one subtraction in the Lippmann-Schwinger 
equation. One should note that for one subtraction n = 1, the dependence on the cutoff is linear. Even in this simple 
example, the renormalization group evolution of the recursive driving term is not straightforward. But it is important 
to stress that one more parameter is required to renormalize Eq. (12 1. Therefore, in agreement with Refs. [571 140) . 
to renormalize the P— wave contact interaction we need two inputs: the scattering volume and the next term in the 
generalized effective range expansion for the cot 6. 

Within our method, with a scaling momentum parameter fi, instead of the cutoff A, the ultraviolet divergence in 
the scattering equation from the potential V{p,p') = \i p p' requires two subtractions (n = 2) to render finite the 
integral appearing in Eq. (10 1. The corresponding result, where we can take the limit A — > cx), is given by 



T{p,p';k^) 
pp' 



Al(fc2,Ai2 



dq 



qHk^ + ii^y 



(M^ 



(fc2 -q2+ ie) 



(13) 



where the parameter Ai is now finite and can be fitted by the scattering volume. However, a dependence on is 
intrinsic to this method, unless the driving term is evolved through a nonrelativistic Callan-Symanzik equation [^,'10]. 

Just to give an impression on how the RG equation works for the P— wave example, we introduce the renormalized 
driving term Vip,p') = V{p,p'; k"^) ~ Ai(fc2, /i^) p p' in Eq. After projection in the P— wave, we have 



9Ai(fc^A^^) 



= 2[Ai(fc2,M')]'- / dqq 

Jo 



(fc2 



(g2+^2)S 



4 



[Al(fc^A^^)]^ 



Therefore, the solution of the RG equation gives 

1 1 



Al(fc2,M') Ai(fc2,Mg) 



2^^ 



3\ 



(14) 



(15) 



Note that, even if at the reference subtraction energy — /Iq the renormalized strength is chosen to be independent of 
/c^, the evolution introduces a fc^ dependence, and implicitly the need of two-subtractions to render finite the T-matrix 
for the P— wave. If one fixes the reference scale at zero energy as suggested in [34], the renormalized strength is the 
scattering volume assuming independence with k^ at this reference energy. In our case, after solving the RG equation, 
the form of the evolved potential is simply the renormalized Ai from Eq. (151 times p p' . Solving Eq. (13) with the 
renormalized coupling from Eq. (15), we get 



k* cot 8 



1 



Ai(fc2,^2) 



(3fc2+^2)^ 



(16) 



The evolution of Ai according to Eq. ( 15 ) shows that it should be linear function of fc^, which demands two constants. 
Indeed, if we re-write Eq. ( 16 ) taking into account the constant plus the k? term from Ai only two independent 
quantities appear, i.e., a scattering volume and an effective momentum. In fact, the form we obtain for k^ co\,b is 
similar to the one obtained in Ref. 1371. 
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III. RENORMALIZATION OF THE NNLO POTENTIAL WITH n = 4 

Once we have established the recursive procedure to renormahze the NN interaction, we need a potential in momen- 
tum space. For the NNLO chiral potential, we adopt a momentum-space form as explicitly given by Epelbaum (19j. 

Even though our method is powerful enough to renormalize the full two-pion exchange (TPE) potential, we consider 
the version with the spectral representation regularization such that the comparison with results obtained by other 
calculations can be more straight. 

For the sake of completeness and the reader's convenience, we repeat here the analytical expressions for the chiral 
NNLO momentum space potential. The LO interaction is given by the one-pion exchange (OPE) plus an contact 
interaction. The strength of the isospin conserving contact terms and their derivatives depends on the LSJ channel, 
considering also angular momentum mixing. We simplify our notation for the contact terms so that the total angular 
momentum and isospin dependence are not shown explicitly in the expressions for the potential. The corresponding 
S'— wave projected matrix elements of the interaction is given by 

VLo{p,p')^VoPEiP,p') + Xo, (17) 
where the unprojected OPE potential, for q = p~ p', is given by 

VoMp,P') = J^[y-J -1-2 + ' 

At NLO, we have some TPE diagrams plus derivative contact interactions. After partial wave projection, it is given 

by 

Vnlo{p,p) = V^p'^^{p,p') + \^{pp')5l,i5l's 

+ (a2 (p2+p'2) + A3 (p' /)) 5l,o<5l',o 

+ A4 (/(5l,2<5l',o+p''^l'.2<5l,o) , (19) 
where the unprojected NLO TPE potential is 

^TPE (P,P ) = - on. 2f4 T^Tv^ i 4M^(5g^ - Ag^- 1) + q (23g^ - IO3.4 - 1) + 



L{q) (ai ■ q) (5=2 ■ q) - q' ffi ■ ^2 } ■ (20) 



V3847r2/4y (27r)3 \ '^^"^'^ J ' h y -^yA ; 4Af2 + g2 

_ / 3gl 

The term proportional to Ai contributes only in the P— waves, the terms proportional to A2 and A3 appear only 
in the 5— waves and the the term proportional to A4 enters only in the coupled channels with j = 1. The A3 term 
(p^ p'^) actually appears at N'^LO in the Weinberg's power counting, but we promoted it to NLO in order to improve 
the S"— waves fit. It could have been promoted to NNLO, but since the iteration of the NLO A2 term (p^ +p'^) leads 
to p"^ p'^ terms, we included it at NLO so that all the p^ p'^ terms combine at once. This has also been done in Ref. 

m- 

Finally, at NNLO, we have other TPE diagrams, with the corresponding unprojected potential given by 



Vnnlo{p,p') = 14pe^°(p,P') 



3si 



{2n)3 \ 167r/4 J \ 16m{AM^ + q^ 



9\ 



+ 



1287rm/4 
5127rm/4 _ 



(ti • T2 



2 

A 

16m 

+ [AMI + 2g' - g\{AMl + 3g')] {2MI + q^)A{q) 



2Ml(2c,-c,)-q\c,^^) 



{2Ml + q^)A{q) 



4M2 



((Ti ■ q){o2 ■ q) — q^{Si ■ S2) 



327r/, 



((Ti ■ q){S2 ■ q) — q [Si ■ S2) 



4 v"^! ■ "^2) 

32) ■ (p ' X p) {2Ml + q^)A{q) 



{2Ml + q^)A{q) 

1 



8m 



(IQMl + 3g') \ A{q) 



647rm/4 



(Ti ■ T2) t (<?! + a2) • (P ' X p) {AMI + q^)A{q) , 



(21) 



7 



where the loop integrals L{q) and A{q) are given by 
L( 



2M^ 



2q 4M2(A2-q2) 
A(q) = — arctan , w e^fA - 2M^) — arctan ~ ^^^^J . (22) 

In the above expressions, A is a spectral regularization scale defined in |19| for the two-pion exchange, in order to 
improve the convergence of the chiral expansion. As suggested A w four pion mass, we set A — 600 MeV. 



A. Evolving the potentials through the recursive renormalization process 

We start by calculating /i^) from the leading order interaction VlOi by solving the Callan-Symanzik equation 

([8| for n = 1 starting at a negative infinite Ji^ up to a finite . The integral form of RG equation Q for n = 1, after 
partial wave decomposition, is 

V^>(p,p\-H) = Vlo{p,P -,-1^ ) + - dqq — " — —^V' ' {q,p ) , (23) 

which brings the leading order interaction to a scale —^^ from its infinitely large fixed-point —Jp j3J [101 [131 US] ■ This 
generates an interaction for the one subtracted scattering equation, which gives the same observables as the leading 
order interaction when /i ^ ^qcd- Now, we obtain V'^'^\p^p'] k^; — /i^) from V^^\p,p'; — /i^): 



V^-\p,p';k-;-f.-)^V^^\p,p';-f.-)-U^ dqq- "^'"'^^frf 1^ ^ fe^; V) 

IT Jo (n^ + q^y 

At the third step we evolve from l^^^) ^j^d introduce the next-to-leading order terms 



(24) 







V^(''(p,p';fc^ V) = ViVLo(p,p'; V) + V'<''(p,p';fc'; V) • (25) 

At the fourth step, the higher order we consider here, we evolve from V'^^'^ and add the next-to-next-to-leading order 
two-pion exchange: 

- - r dq q^ V^'\p, q- fc^; -n^) + % V('\q,p'- k'; -f,') 

71" Jo (p + r) 

v('\p,p'-e-~^i^) = VMNLo{p.p'■-^^^) + v^^\p,p'■,e■-p^) . (26) 

With the above, the half-on-shcll T— matrix with four subtractions is a solution of 

TT Jq \P + P J K - P + le 



(27) 



Note that in the above equation the term given by [{p^ + k'^)/{p'^ + p'"^)]* works effectively as a regulator, canceling 
the singularity presented in the starting interaction. The driving term is generated by consecutive subtractions in 
the kernel at some defined energy scale fP, with the advantage that such scale can be moved freely as long as V^^^ 
satisfies the RG equation, given by Eq. (Isl), that guarantees that the scattering amplitude is unaltered. 
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It is important to explain the strategy used to insert the potential in the recursive to obtain the driving term V^^^ 



that enters in Eq. (27). We include each order of the interaction (LO, NLO, NNLO) in the step where there is enough 
subtractions to renormalize it. Consequently, we insert the LO potential in the first subtraction, the NLO potential 
in the third subtraction and so on. Tables |T] to |IV| list the parameters for the potentials used in this work. 

Before closing this subsection, we discuss briefly our scheme and the one used in Ref. ^34J, where S*— waves have 
been considered with more than one subtraction. In that work, the successive subtractions are performed in a rather 
different form compared to our approach. They call it a mixed scheme. The main difference is: in order to handle the 
momentum-dependent contact interaction + p'^, the cutoff also plays an important role in |34j . since according to 
our method this interaction actually requires three subtractions. In our case, we can approach any kind of interaction 
as long as we perform enough subtractions. The cutoff has no physical relevance since the multiple subtractions allow 
us to sum up to arbitrarily large momentum values. 

1 c 

TABLE I: Strengths of the LO contact interactions, which reproduce the scattering lengths for the S waves. The values of Aq " 
and Ag'^^, in units of fm, are given at the energy scale —p,^, with p, = 30 fm^^ (/i^ = 41.47x 900 MeV). 



Strengths 






Ao(fm) 


-0.0203 


-0.24142 



TABLE IL Strengths of the contact interactions for the fits with the LO potential plus the NLO contact interactions. The 
values of Ag'^" and Ag'^^ are given at the same energy scale as in Table|l](-^i2 = -41.47X 900 MeV); with Aa^" and x'/° at 
-^i^ = -50 MeV. The other strengths are given at -jU^ =-100 MeV. 



Strengths 


'So 'Po 


'Si 'Pi 'Pi 'P2 ei 


Ao(fm) 


-0.0165 - 


-0.2480 - - - - 


Ai(fm3) 


0.25 


0.04 0.007 -0.07 - 


A2(fm3) 


2.2660 - 


0.1 - - - - 


A3(fm'^) 


2.0047 - 




A4(fm=') 




- 0.001 



TABLE III: Strengths of the contact interactions for the fits with the full NLO potential. The values of the A's are given for 
fL^ and /i^ as in Table [HI 



Strengths 


'So 'Po 


'Si ^Pi 


'Pi 'P2 ei 


Ao 


-0.0190 - 


-0.1602 - 




Ai(fm^) 


0.37 


0.063 


-0.078 -0.04 - 


\2{im') 


2.2660 - 


0.1 




A3(fm5) 


2.0047 - 






A4(fm^) 






- 0.17 



TABLE IV: Strengths of the contact interactions for the fits with the NNLO potential. The values of the A's are given for 
and fj^ as in Table [ill 



Strengths 


'So 'Po 


'Si 'Pi 


'Pi 'P2 tl 


Ao(fm) 


-0.0189 - 


-0.1217 - 




Ai(fm=') 


0.303 


0.066 


-0.19 -0.1 - 


A2(fm3) 


2.2660 - 


0.1 




A3(fm5) 


2.0047 - 






A4(fm') 






- 0.17 
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FIG. 1: (color online) Half-on-shell recursive driving terms 1/'"' (g, fc; fc^) for the ^So channel, corresponding to the energies 
E = {hk)'^ /mN = 50, 100, 150 and 200 MeV (with units such that h — 1 and mjv — 1). The lines identification, shown inside 
the lower-left panel, are given for all four panels. The values of the strengths, as well as the subtraction energy scales (not 
explicitly shown here), are given in Tables [I|-[IVl 



B. Half-on-shell amplitudes for the recursive renormalization process 



Now, let us consider the calculation of the half-on shell matrices F*^"^ {q. k) for n =1 to 4, from the solution of ( 23 ) to 
(26), and the half-on-shell matrix elements of r(g, fc; fc^) solution of (27 1. Before starting our discussion, a side remark 



is worthwhile. As we have seen before, from the above equations, the recursive driving terms obey integral equations 
with kernels defined by subtracted free Green's functions multiplied by recursive driving terms one order below. From 
the subtracted method itself the lowest order recursive driving term satisfies an equation which is energy independent. 
In the following orders the integral equations for the driving term depend on fc^ , and the subtracted Green's function 
is negative, which brings a curious effect on the respective solutions of ([6|. The corresponding homogeneous equation, 
in the case of attractive recursive interaction, suffers an enhancement by increasing fc^, due to the factor (/x^ -I- k"^) 
inside the kernel, which may allow it to have a solution. When this happens, an unphysical pole will occur in the 
recursive interaction of the integral equation. This can be realized even with a regular attractive potential if multiple 
subtractions are used to compute the scattering amplitude. This unphysical pole is completely washed out in the 
solution of the n— subtracted LS equation given below. 

We can observe the half-on-shell recursive driving terms evolving through the four subtractions from V^^\q^ k) up 
to t{q,k) for some values of the on-shell momentum k in Figs. [l]to [4j for the ^Sq and '^S'l channels. In Figs, [l] and 
[3] we have the driving terms and half-shell T-matrix elements, respectively, for the ^5*0 channel. Correspondingly, 
we have the results for the ^Si channel in Figs. [2] and |4j These figures exhibit us an interesting finding due to 
Redish and Stickbauer [JT]. They observe that distinct half-on-shell potentials may lead to similar half-on-shell 
T— matrix. In other words, given different V{q, k) which fit the same on-shell observables their corresponding half- 
on-shell T— matrix should be quite equivalent despite the V{q, k) discrepancies. That is why we observe a smooth 
behavior of the scattering amplitude with energy in Fig. [2] while the recursive driving terms, in particular V^^\ vary 
considerably as can be seen in the right panel of Fig. ([T]). 



The set of integral equations for the subtracted driving terms Eqs. (23)- (26 1 given above also deserves further 



comment to remove the naive misconception that the approach given in [10 



"involves invoking the Born approximation 



and consequently unreliable for the higher singular potentials" 34]. As one should observe, such approach relies on 
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the definition of a driving term F^") at each, subtraction order. Of course, such term is the transition matrix at the 
subtraction point. The driving term should be evolved dynamically by the renormalization procedure to the next 
order, where higher divergent interactions are added, if necessary. If a contact interaction is included, its strength is 
the renormalized one at that point. Such procedure continues up to the number of subtractions required to render 
finite the corresponding scattering equation. Let us emphasize that, differently from the usual Born approximation 
(where the dynamics does not evolve the potential), this renormalization procedure relies on a recursive evolution of 
the driving term, from one order to the next one. Moreover, it produces a unitary S"— matrix. 

The LO term, OPEP plus Dirac-(5 for •^5'i Di and and OPEP for the higher partial waves are brought from 
the fixed point at infinite to the reference scale where the physical information is supplied to the two-nucleon 
system. Moreover through the renormalization group equations we are able to evolve precisely the driving term of 
the subtracted equations to any arbitrary scale without altering the physical content of the observables. However, in 
this case the operator form of will acquire a non-trivial form and will be not easily identified with the starting 
potential. 



C. Full-on-shell amplitudes along the recursive renormalization process 



The next step is the computation of the full-on-shell matrices ^'■"•'(fc, k; k^; — /i^), obtained from the solution of the 



recursive Eqs. (23) to (26), with p = p' = k. Once we have the on-shell T— matrix, computed from Eq. (27) with 



p = k, we can obtain the S*— matrix as 

S = l-2ikT. (28) 

Figures [5] -[7] present the full-on-shell recursive potentials, as we step forward in the recursive subtractive renormal- 
ization process, and the real part of T(fc, k; k"^) (denoted by T(fc, k) in the figures) in the uncoupled states with total 
angular momentum up to j = 2. For coupled channels, with j — 1 and j = 2, the evolution of the recursive potential 
and the real part of the scattering amplitude are shown in Figs. [8] and |9] 

Results for ^Sq and ^Pq are shown in Fig. [H] Clearly the results of scattering amplitude for the singlet S*— wave 
channel are far from perturbative. At zero energy, it is observed a strong deviation of r(0, 0; 0) (that gives the 
scattering length) from ^(''^(0, 0; 0; — ^^), which is due to the nearby singlet virtual state. The contribution of the 
contact in NLO to ^(■^^(fc, fc; A:^; — /i^) distinguishes it from the recursive process. In particular, ^(''^(fc, fc; /c^; — /i^) 
presents a pole, which arises from the solution of the integral equation, that does not affect the scattering amplitude 
in the energy range that is shown. The comparison of V^'^'>{k,k;k^;—fi^) with the real part of T{k,k;k^) for ^Pq 
shows a behavior that indicates the dominance of the driven term. The zero of both quantities are close together and 
the magnitude of the real part of T(fc, fc; k^) decreases in respect to ^^'^•'(fc, k; fc^; — /i^) with energy. The contribution 
of the NLO contact to V^'^''{k, fc; k^; — /i^) is also clearly seen in the figure, while the TPE NNLO potential appears to 
be not so much relevant in this wave. 

The on-shell matrix elements of the recursive potential for ^Pi and '^Pi and scattering amplitude are shown in 
Fig. [6) For these waves the TPE NNLO potential does not gives a relevant contribution as one observe from the 




q [fm '] 



FIG. 2; (color online) Half-on-shell T-matrix elements for the ^So channel, for the energies = 50, 100 and 200 MeV. As in 
Fig. [Tj the values of the strengths, as well as the subtraction energy scales (not explicitly shown here), are given in Tables |l]- 
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FIG. 3: (color online) Half-on-shell recursive driving terms 1/'"' (g, A:; fc^) for the "^Si channel, corresponding to the energies 
= 50, 100, 150 and 200 MeV, as given inside the upper-right panel. As in Fig. [l] the values of the strengths, as well as the 
subtraction energy scales (not explicitly shown here), are given in Tables |I]-|IV] 



difference between V'^^\k, k; k^; — /i^) and V^-^^{k, k; k"^; — A*^)- The NLO contact is important in these waves, as also 
seen for ^Pq. The real part of scattering amplitude for ^Pi is also somewhat dominated by V^'^\k,k;k'^;—iJ?). The 
uncoupled D— waves are shown in Fig. [7] No contacts are present and the recursive potentials change smoothly from 
one to the next. 

The coupled ■^S'l Di on-shell potentials and scattering amplitudes are presented in Fig. [s] The real part of the 
scattering amplitudes for these channels are far from perturbative. At zero energy the real part of t{0, 0; 0) for the 
■^S*! channel, that gives the scattering length, moves strongly from V'^^^ {0,0;0; —^j?) due to the deuteron pole. The 
same is observed at low energies for the off-diagonal amplitude, related to the importance of the deuteron D/S ratio 
to the mixing parameter [32]. The contribution of the contact in NLO to V''^^{k, k; k^; —fJ,^) distinguishes it from the 




FIG. 4: (color online) Half-on-shell T-matrix elements for the ^Si channel, for the energies fc^ = 50, 100, 150 and 200 MeV. As 
in Fig.rn the values of the strengths, as well as the subtraction energy scales (not explicitly shown here), are given in Tables m 
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E = k" (MeV) 

FIG. 5: (color online) On-shell elements V^"'(fc, k; fc^) and T{k, k; k^) for the uncoupled channels with J = 0, as function of fc^. 
The legends for the curves, given in the upper panel, are the same for both panels. 



recursive process. The pole of V^-^\k, fc; k^; — M^) does not affect the scattering ampUtude in the energy range of our 
calculations. The coupled — F2 on-shell potentials and scattering amplitudes are seen in Fig. and we observe 
a strong increase in the magnitude of V'^^^k, k; fc^; — /i^), that also is an indication of the contribution of the NNLO 
potential to this wave. 





100 125 
E = (MeV) 



FIG. 6; (color online) On-shell elements V'"'(fc, k\ k^) and T{k, k; k^) for the uncoupled channels with J : 
The legends for the curves, given in the upper panel, are the same for both panels. 



1, as function of k^ . 
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FIG. 7: (color online) On-shell elements k; k^) and T{k, k; k^) for the uncoupled channels with J = 2, as function of k^. 

The legends for the curves, given in the upper panel, are the same for both panels. 



IV. RESULTS FOR NUCLEON-NUCLEON PHASE-SHIFTS AND MIXING PARAMETERS 



For the analysis of the phase-shifts and mixing parameters obtained with the renormalized strengths and subtraction 
energies presented in Tables [l] to IV we adopt a systematics which splits the calculations in four sets: 

(i) leading order (LO), as given in Eq. (17); 

(ii) leading order plus next-to-leading order contact interactions (LO -I- NLO CI), as given in Eq.(19); 

(Hi) full next-to-leading order (NLO), consisting of one-pion exchange, two-pion exchange at NLO, and contact 
interactions; 

(iv) next-to-next-to leading order (NNLO), which is the NLO plus TPE diagrams at NNLO. 

In particular, the set (ii) was inspired in the idea to promote some NLO terms to LO, as a way to overcome difficulties 
with the Weinberg power-counting rule (see Ref. [3]). One consequence of the failure of the naive dimensional analysis 
(NDA) is that contact interactions that are subleading in Weinberg's power conting are in fact necessary in order to 
make sense the T matrix even at LO. In Refs. [2H1|, it is shown how the origin of the failure of NDA is the singularity 
of one pion exchange (OPE), for which a single counterterm suffices at LO. 

For the LO potential in the singlet and triplet channels, the respective scattering lengths are fitted. Indeed this 
calculation reproduces the results obtained in Ref. ^ , which we supply here for completeness and in order to compare 
with the results obtained in NLO and NNLO. The results for the phase shifts for the waves ^Sq and '^Pg are shown in 
Fig. 10 We should observe that results corresponding to the LO -I- NLO CI were already presented in Ref. [lU]. The 
NLO and NNLO calculations were done by using the same renormalized strengths for the derivative of the contacts 
(see A2 and A3 in Tables [ll] - IV) and only refitted the singlet scattering length by changing slightly Aq. This means 
that at the range of energies we perform our calculations the contributions of TPE potentials in this wave are small. 
Moreover, the NLO and NNLO calculations present a small systematic deviation steadily increasing with energy. 
This is possibly due to the strong attraction of the corresponding TPE NLO and NNLO potentials which increases 
at higher momentum. 

In the "^Pq shown in the lower frame of Fig. 10 there is a contact interaction of the type Vcontact — ^1 P p' at NLO 
along with one-pion and two-pion exchanges. The LO calculation is not able to fit the scattering volume of "^Pqi while 
with Ai given in Tables |ll] to IV we obtain a better fit as shown by the corresponding adjustment of the low-energy 
Nijmegen phase-shifts. It is seen that a small change in Ai is enough to give the scattering volume for sets (ii), (Hi) 
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FIG. 8: (color online) On-shell elements 1/'"'(fe, k; k^) and T{k, k; k^) for the coupled channels with J = 1. In the middle panel, 
the amplitude T(fc, k\ k^) for t\ has been multiplied by a factor of 10 in order to highlight its qualitative behavior. The legends 
for the curves appearing in the three panels are the same as the ones given in Fig. [7] 



and (iv). The small change in the contact term which fits the singlet scattering length is similar to what is found for 
^Pq in respect to the scattering volume. 

The phase shifts for the waves ^Pi and ^Pi are shown in Fig. 11 In both channels there is a contact interaction of 
the type Vcontact — P p' at NLO along with one-pion and two-pion exchanges. In this case, we observe that the LO 
+ NLO CI is already giving a quite good fit, when comparing with the Nijmegen results. The TPEP contribution is 
marginal as observed by the slight change of Ai required to keep the fit. (See Tables [lT| to |IV|) . 

The results for the phase shifts, for the uncoupled channels and ^-D2, are shown in Fig. 12 As already mentioned 
before, no contact terms are presented, in this case. 

In the following, we discuss the results for the coupled channels for the spin triplet j — 1 in the '^Si Di states 
and j = 2 in the ^P2 F2 states. The phase shifts for the ^Si Di coupled channels are shown in Fig. 13 and for 
F2 in Fig. [l4| 

For the mixing parameter ^Si Di, we consider ei, as defined in Ref. [13], instead of the Blatt-Biedenharn 
definition, cbb IS]- The reason for this choice is related to the precise measurements available for ei. We observe in 
Fig. 13 that the mixing parameter ei can be well fitted with only the leading order plus a small contact term in the 
mixed states (see A4 in Table |ll|. This indicates that the physics of ei seems well controlled by OPEP, as long time 
ago predicted by a shape-independent expansion [32] • 

What do we see when the NLO and NNLO potentials is inserted in our method? The aforementioned nice fitting 
of ci by the LO potential plus small contact term in the mixed channel disappears [see Fig. 13 . We met the widely 
recognized difficulty that the effective potential has problems in the describing ei. To make concrete this point we 
mention that different renormalization approaches, coordinate space renormalization method |30| and subtraction 
plus cutoff [31] also faced the same difficulty to fit ei in NLO and NNLO. In particular Ref. [31] exploited the strong 
momentum cutoff dependence to fit ei, however the fit is not robust in the sense that it should be smoothly cutoff 
dependent. This common difficulties in different and independent calculations reveals the distinct role played by the 
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FIG. 9: (color online) On-shell V^"\k, k; k'^) and T{k, k; k^) for the coupled channels with J = 2. The legends for the curves 
appearing in the three panels are the same as the ones given in Fig. [7] 

singularities in different waves. 

From the results shown in Fig. [TSj we see that the TPEP does not exhibit a systematic behavior in the different 
partial waves regarding the NLO and NNLO potentials. From our point of view a systematic behavior would require 
a cutoff for TPEP or its inclusion only after higher order singular terms are included. This problem is acute in the 
mixing parameter of the coupled channel ^Si Di. The mixing parameter with OPEP plus NLO contact is well 
fitted up to 100 MeV lab energy. The introduction of NLO and NNLO TPEP, clearly worsen the fit to the Nijmegen 
phase-shift analysis |45| . To fit the mixing parameter it is important to have the deuteron asymptotic D/S ratio within 
their accepted value. This was pointed out in a recent work |31j . where the authors have also included the contribution 
of the Delta excitation. Indeed, long-time ago, it was already concluded that the mixing parameter at low energies 
is determined by the deuteron properties and by OPEP By considering the Blatt-Biedenharn definition, 

it was demonstrated in that the correct long-range behavior of the tensor potential is essential for a realistic 
reproduction of the mixing parameter. As shown in j42j . separable tensor Yamaguchi and square- well potentials, 
which do not possess the OPE tail, when fitted to reproduce the deuteron binding and asymptotic normalization, 
badly fail to reproduce the correct value of the mixing parameter. The NLO and NNLO potentials for the deuteron 
channel seem to give a too strong contribution, enhanced by the singular behavior of OPEP at short distances in this 
channel. Although, the diagonal and off-diagonal NLO contact plus OPEP and LO contact are enough to give a nice 
fitting to the mixing parameter up to E'lab. — 100 MeV, the inclusion of NLO and NNLO TPEP does not provide 
a good fit, and the results systematically deviate from the Nijmegen data for ei. This indicates that by keeping the 
OPEP potential intact, the NLO and NNLO TPEP potentials have to have a cut at short distances. 

It is appropriate to summarize the comparison between our method of multiple subtractions and that used in ref. 
[34j . In ours no cutoff is needed while in the later one just one subtraction is required since higher order singularities 
coming from the delta derivatives are tamed by a cutoff. Numerically, in general, both procedures lead to a similar 
fitting for the S and P waves up to energies of about 200 MeV. We also remark that the mixing parameter, ei, presents 
the same deficiency in the fit by both methods as well as for a different regularization and renormalization method in 
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FIG. 10: (color online) Uncoupled channels for J = 0. Phase shifts for the ^So wave with the subtraction point at -50 MeV, 
and for the ^Pq wave with the subtraction point at -100 MeV. The legends for the curves, given in the upper panel, are the 
same for both panels. 

coordinate space without the contribution of the Delta excitation [3T] . 

We note that, once the parameters are fitted (renormalized strengths and subtraction points) the subtraction 
scale can be moved without changing the T-matrix results. Our method is RG invariant by construction. The flow 
equation ([s]), which transforms the driven term of the subtracted scattering equation (27 1 under the dislocation of 
the subtraction point, was derived from the invariance of the scattering matrix under variations of \x. 

To conclude this section, we illustrate the cutoff independence of our method. We choose to present results for two 
singular cases, the ^5*0 and ^Pq phase-shifts at 50 MeV. They were obtained by considering the NNLO potential, with 
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FIG. 12: (color online) Phase shifts for the ^L»2 and ^L»2 waves. The legends for the curves are given in the upper frame for 
both the cases. 



4 subtracted scattering equation (27). As shown in Fig. 15 the phase-shifts stabihze when the cutoffs are driven 



towards infinity. The same kind of cutoff behavior is verified for the other orders of the potential (LO and NLO), 
as well as for other partial waves and Eiab energies. Our results shown along this paper were obtained for infinite 
cutoffs, with the verified convergence due to the subtractive renormalization approach we are considering. 



CONCLUSION 



We present a systematic methodology to renormalizc the nucleon-nucleon interaction using a recursive subtracted 
approach with multiple subtractions in the kernel. Within the subtractive scheme we studied the two-nucleon scat- 
tering T— matrix for the NLO and NNLO potentials, for all partial waves up to j = 2. The renormalized strengths 
of the contact interactions, the so-called low-energy constants, were fitted to the low-energy phase-shifts and mixing 
parameters from the Nijmegen partial wave analysis for a reference subtraction point. 

In order to show how the multiple subtraction renormalization scheme can be implemented, an analytical example 
is given for the P— wave channel. In that case, two subtraction are required to eliminate the ultraviolet singularity of 
the interaction. We give explicitly the solution of the RG equation. Although we have not fully explored the analytical 
form of this P— wave amplitude, we call the attention of the reader to the richness of the analytical continuation to 
the complex energy plane in order to obtain virtual and resonant states [371 SHI HZ] • 

In a sharp or smooth cutoff approach, the LO, NLO and NNLO potentials are regularized to vanish above a certain 
momentum scale and then inserted in the LS equation. In our method, the cutoff is purely instrumental and the limit 
of the momentum cutoff going to infinity can be easily performed, since a finite T— matrix is ensured by multiple 
subtractions in the scattering equation. The original potential is kept intact and enters in the recursive process as we 
described in detail. 

We analyze the matrix elements of the potentials in momentum space for the ^Sq channel as we go from LO through 
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NLO up to NNLO. This suggest in a practical way a momentum scale, which we associate with the subtraction point, 
where the systematic expansion of the potential from the chiral effective field theory should be used within our 
subtraction scheme. From such analysis, we verify that a momentum scale of about 1 fm^^ separates the matrix 
elements of NLO and NNLO potentials in two regions: below this scale they are comparable to the OPEP. With 
this indication, we choose a subtraction point at an energy — /x^ ~ -50 to -100 MeV. The subtraction momentum 
scale comes to be at the order of the QCD scale, Aq(;7£)[35] and well below the p— meson mass, consistent with the 
general suggestion of Weinberg [T]. Important to observe that, by taking into account the values of the NN low energy 
parameters and the form of the bare interaction including the derivatives of the contacts, a renormalization scale 
near the value of the h^qcD is consistent with the Wigner bound [321 ISO] ■ Once we fit the renormalized strengths of 
the contacts, considering the pion exchange potentials at LO, NLO and NNLO, the subtraction point can be moved 
arbitrarily. Our method is RG invariant by construction, and the fiow equation, ([8|, transforms the driven term of 
the subtracted scattering equation such that scattering matrix invariant under dislocation of /i. 

We show how the half-on-shell potentials for ^Sq and ^Si channels evolve through the four subtractions, from 
fc) up to T{q,k;k^). These exhibit a relation with the interesting finding of Redish-Stickbauer [5T] that the 
half-on-shell potential and corresponding half-on-shell T— matrix can be very different. Given distinct V{q,k), which 
fit the same on-shell observables, the corresponding half-on-shell T— matrix should be quite equivalent in spite the 
sharp V{q, k) differences. That is consistent with the smooth behavior of the scattering amplitude with energy, while 
the recursive driving terms can vary considerably. 

In the ^Po channel, a contact interaction is introduced at NLO together with OPEP and TPEP. The LO calculation 
is not able to fit the scattering volume of ^Pq, while the addition of contact gives a better fit of the Nijmegen phase- 
shifts. It is verified that a small change in the renormalized strength of the contact, obtained only with OPEP, is 
enough to reproduce the scattering volume for NLO and NNLO. The same behavior is observed for the strength of 
the contact interaction when the singlet scattering length is fitted with NLO and NNLO potentials. 

For the ^Pi and ^Pi channels we introduce a contact interaction at NLO along with OPEP and TPEP. A quite 
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good fit of the Nijmegen phase shifts is obtained with OPEP plus the contact. Again we observe that the TPEP 
contribution is marginal as a slight change of the renormalized strengths keeps the fit. These observations indicate the 
absence of an essential singularity in the OPEP and TPEP potentials up to NNLO in uncoupled P— waves, beyond 
the contacts itself, differently from what is found for the coupled ^Si Di channel. 

After our fitting of the contact interactions for the P— waves (^Pq, ^Pj, ^Pi, ^-Pi), the following observations apply 
to the matrix elements in momentum space of the LO, NLO and NNLO potentials. For "^Pq and ^Pi, the contact 
dominates above 3 fm~^, while below 2 fm~^ the LO potential dominates. For the '^P2, the LO potential is weak, 
while the inclusion of the contact interaction enhances the attraction of the NLO potential, which weakens the strong 
repulsion from the NNLO TPEP. 

In the next, we summarize our findings for the coupled channels: ^ Si Di and '^P2 F2. The mixing parameter 
ei can be well fitted at LO plus a small contact term which couples the S — D states in the interaction. This indicates 
that the physics of ci seems well controlled by OPEP, as suggested in [42" by a shape-independent expansion. When 
we include the NLO and NNLO potentials within our method, the nice fitting ei is destroyed. The NLO and NNLO 
potentials for the deuteron channel seems to give a too strong contribution, enhanced by the singular behaviour of 
the OPEP at short distances in this channel. 

The difficulty in fitting ei with effective potentials in NLO and NNLO, is recognized by different renormalization 
approaches [301 IMj- Such shortcoming is related to the strong momentum cutoff dependence of ei, as found for 
example by the coordinate renormalization approach [30] , possibly due to the singular behavior of OPEP [3S] . 

The contribution of TPEP does not exhibit a systematic behavior in the different partial waves regarding the 
NLO and NNLO potentials. Our results suggest that, for a systematic behavior, the TPEP should be weakened for 
momentum larger than few fm~^'s. This could be done via a cutoff or by considering higher order contacts, which 
could suppress the TPEP contribution in the appropriate partial waves. 

We found that the derivative contact terms dominate over the NLO and NNLO two-pion exchange interactions 
in the S*— wave channels, starting at low-momentum scales (~ 0.3 fm~^). In P— wave channels the contacts are also 
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FIG. 15: Phase shifts for the ^Sq (left panel) and ^Pq (right panel) waves, for Eiat ~ 50 MeV, exemplifying a general behavior 
of the phase shifts for increasing cutoffs. Both the cases presented here were calculated considering the NNLO potential, using 
the n =4 subtracted scattering equation (27 1, with parameters given in Table llVl The same kind of behavior is verified for the 



other orders of the potential (LO and NLO), as well as for other waves and Eiab energies. 



important for the fitting of the corresponding phase-shifts. 

Finally, we should observe that the input of our method is the T— matrix at a given energy, where the physical 
information is supplied to the two-nucleon system. When the interaction is at LO, the point where the physical input 
is given is not constrained. Once we move to NLO and NNLO, the energy — /i^ arises as a scale where the low-energy 
observables can be obtained. In view of that, the particular value of the subtraction point acquires the status of a 
physical scale when the NLO and NNLO interactions are introduced. It is gratifying to verify that our fittings with the 
associated subtraction point, given by the renormalization scale, comes to be about Aqcd, well below the p— meson 
mass. This is consistent with the general Weinberg's[T] suggestion, that an effective potential should be valid for 
momenta much smaller than a typical QCD scale of 1 GeV, and the intermediate nucleon-nucleon propagation should 
be damped at such small momentum scale. 
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Appendix A: NN Phase Shifts and mixing parameters 

We follow the definitions given in Ref. [43 for the nucleon-nucleon phase shifts and mixing parameters, which are 
appropriate to the case that we have coupled channels. 

Considering the cases with j > 0, for the coupled channels we have an S*— matrix expression given by 

_ / \( e'^-i \ / cos(2ej) isin(2e,) \ ( e'^-^ \ 

from where the corresponding phase shifts and mixing parameter are 
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Oj + l 



2 \Re ij+i j+i , 



(A3) 



£„■ = — tan 
^ 2 



^ ('S'i-ij+i + ■S'j+i.j-i) 



(A4) 
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